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ABSTRACT 


Objective: The monthly mean precipitation over India is used to investigate natural periodic wave characteristics using a novel 
technique of Empirical Mode Decomposition (EMD) and Ensemble Empirical Mode Decomposition (EEMD). 

Methods: These two methods empirical mode decomposition and ensemble empirical mode decomposition are interesting ap¬ 
proach to decompose the signals into locally periodic oscillations. 

Result: The Intrinsic mode functions (IMFs), will easily Identify the Embedded structures, even with those smaller amplitudes. 
Ensemble Empirical Mode Decomposition is better performed than the Empirical Mode Decomposition technique. The Empirical 
Mode Decomposition method was observe the mode mixing of two signals in Intrinsic Mode Function, but Ensemble Empirical 
Mode Decomposition found the distinct and clear peak of each periodic signal in each Intrinsic Mode Function. 

Conclusion: In Empirical Mode Decomposition method has 11 Intrinsic Mode Functions and Ensemble Empirical Mode Decom¬ 
position has 10 Intrinsic Mode Functions, it is due to the Noise-assisted method to reduce the Intrinsic Mode Function numbers. 
The decomposed oscillations in Ensemble Empirical Mode Decomposition are above confidence interval and significant, it has 
515 iterations than Empirical Mode Decomposition method has 740 iterations. We observe the computational time is lesser in 
Ensemble Empirical Mode Decomposition than Empirical Mode Decomposition method. 

Keywords: Precipitation, Empirical mode decomposition (EMD), Ensemble empirical mode decomposition (EEMD), Periodic oscil¬ 
lations, Lomb-Scargle (LS) spectral analysis 


INTRODUCTION 

Precipitation is probably the most important of the essential 
climate condition and its crucial role to sustain any form of 
life on earth as a major source of fresh water, its major im¬ 
pact on weather, climate change, and related issues of soci¬ 
ety’s adaptation. The occurrence of precipitation is highly 
variable in space and time. Finally, high-quality monthly 
precipitation data sets across a long-term period are key in¬ 
formation for an improved understanding of the global water 
cycle (Becker et al., 2012). 


The spatial and temporal variations of rainfall are important 
in understanding the hydrological balance on regional and 
global scales. The distribution of precipitation is also criti¬ 
cal for water control in agriculture, power generation and 
drought-monitoring. Nishant Malik et al., (2011) evalu¬ 
ates the Indian summer monsoon (ISM) rainfall over South 
Asia is the result of the interaction of several complex at¬ 
mospheric processes evolving at many different spatial and 
temporal scales (e.g., Webster 1987). By the influences of 
the interplay of synoptic scale weather phenomena, the In¬ 
dian summer monsoon rainfall patterns are also modulated 
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by the steep topography of the Himalayas (e.g., Bookhagen 
and Burbank 2010). Hence, monsoonal rainfall has highly 
complex spatiotemporal patterns. 

The Indian summer monsoon (June to September) rainfall is 
essential for the economic development of population, dis¬ 
aster management, hydrological planning for the country by 
Guhathakurta and Rajeevan (2008). Earlier, Parthasarathy 
(1993 and 1994) used 306 uniformly distributed rain-gauge 
stations for construct the precipitation series. Attempts have 
been to study the annual, seasonal and long-term trends for 
the Indian region as well as for smaller sub-divisions using 
rainfall data (Parthasarathy et al., 1993 and P. Kishore et al. 
2015). 

It is well known that, the rainfall during monsoon season 
over Indian region exhibits, large spatial, temporal, intra- 
seasonal and inter-annual variability. It is evident that Indian 
Summer Monsoon (ISM) exhibits different variations with 
different periodicities starting from active and breaks pe¬ 
riod too, intra-seasonal, inter-annual, quasi-biennial oscilla¬ 
tion (Rao and Lakhole, 1978), El~Nino Southern oscillation 
(Shukla and Paolino, 1983), solar cycle (Bhalme and Jadhav, 
1984). Their analysis reveals that the annually sampled sea¬ 
sonal data is characterized by near periodic oscillations of 3, 
5.8, 11.6, 20.8 and 37 year periods. In general, this variabil¬ 
ity is extracted from long-term data sets by using different 
methods like Fourier and wavelet analysis. From the above 
analysis, the authors have concentrated the only one oscilla¬ 
tion but not all. 

In this study, we make use of two different methods like em¬ 
pirical mode decomposition (EMD) and ensemble empirical 
mode decomposition (EEMD) to extract different oscilla¬ 
tions present in the long-term rainfall data over Indian re¬ 
gion. 

Empirical mode decomposition is a form of adaptive time 
series decomposition method by Huang et ah, (1998). Some 
standard forms of spectral analysis methods like Fourier 
analysis assume that a time series (either linear or nonlinear) 
can be decomposed into a set of linear components. 

In contrast, the Empirical mode decomposition method does 
not assume a time series is linear or stationary before anal¬ 
ysis, it lets the data speak for themselves. Empirical mode 
decomposition adaptively decomposes a signal into a set of 
intrinsic mode functions and a residual component. When 
the intrinsic mode functions and residual are summed to¬ 
gether, they form the original time series (Srikanthan et ah, 
2011). An inconvenient feature of Empirical mode decom¬ 
position is mode mixing, where a fluctuation of given fre¬ 
quency may split across two intrinsic mode functions (Peel 
et ah, 2011a). The adaptive iterative nature of the Empirical 
mode decomposition algorithm means mode mixing is dif¬ 
ficult to avoid without subjectively deciding on the likely 



nature of any signal to be extracted before analysis. Mode 
mixing between intrinsic mode functions is problematic, to 
investigate the significance of intrinsic mode functions, as 
an expected physical signal may be present but split across 
intrinsic mode functions. Wu and Huang (2009) proposed as 
Ensemble Empirical mode decomposition, it is a noise-as¬ 
sisted data analysis method, to overcoming the mode mixing 
problem in intrinsic mode functions. 

In Ensemble Empirical mode decomposition, an ensemble 
of Empirical mode decomposition trials is obtained by add¬ 
ing white noise to the time series of before the each Empiri¬ 
cal mode decomposition run. The intrinsic mode functions 
and residual from each trial are grouped by intrinsic mode 
functions order into ensembles, and the intrinsic mode func¬ 
tion and residual ensemble averages to form the Ensemble 
Empirical mode decomposition. Since the white noise is dif¬ 
ferent for each trial of Empirical mode decomposition and 
its noise cancels out during averaging as the ensemble size 
increases. However, the noise serves the useful purpose of 
changing the order of local maxima and minima within the 
time series, thus different Empirical mode decomposition 
outcome in each trial is formed. Wu and Huang (2009) be¬ 
lieve Ensemble Empirical mode decomposition method pro¬ 
vides more physically meaningful intrinsic mode functions 
and residue than the traditional Empirical mode decomposi¬ 
tion method. 

We present here different period of oscillations using the 
Indian Meteorological Department precipitation data. The 
data from 1901-2010 has been used in the present study. The 
analysis is expected to provide more information in the peri¬ 
odic oscillations using different spectral analysis techniques. 
In the next section, a brief description of our data collection 
of analysis procedure is given below section2. Results and 
discussion are given in the subsequent section3. Finally, our 
results are summarised in section 4. 


DATA AND METHODOLOGY 
IMD Precipitation Datasets 

The Indian Meteorological Department (IMD) l°xl°gridded 
precipitation datasets of the periods from 1901-2010 over the 
Indian region (Rajeevan et al., 2006, 2008) is used for the 
present study. This data analysis and results are organized 
from 3700 rain-gauge stations over India. Each grid consists 
of several stations of data and linear interpolation technique 
(Shepard, 1968) is used to provide the missing data points. 

Empirical Mode Decomposition (EMD) tech¬ 
nique 

The Empirical mode decomposition method was firstly intro¬ 
duced by Huang et al., (1998). The essence of the approach 
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is to empirically perceive the intrinsic oscillatory modes by 
way of their functional time scales within the statistics to de¬ 
compose them consequently. Empirical mode decomposition 
method makes use of local features time scale of the signal, 
extracting some intrinsic mode functions and residual from 
the original signal and intrinsic mode function show the local 
features of the data while the residual component shows the 
slow change of the signal. The key idea of this method is em¬ 
pirical mode decomposition, and it can make any complex 
data sets be decomposed for a limited, usually a few num¬ 
bers of intrinsic mode functions. An intrinsic mode function 
meet two conditions: (1) In the whole data set, the number of 
extrema and the number of zero crossings must either equal 
or differ at most by one. (2) At any point, the mean value of 
the envelope defined by the local maxima and the envelope 
defined by the local minima is zero. It is versatile in a broad 
range of applications for extracting signals from data gener¬ 
ated from non-stationary processes (see an example, Huang 
et al., 2008; Kishore et al., 2012) the decomposition proce¬ 
dure is as follows. 

Let the original time series of monthly precipitation data be 
P o (t). First the upper (E+) and lower (E-) envelopes of local 
maxima and minima respectively are estimated using a cubic 
spline interpolation. Next, at each time instance, the mean 
value of the two envelopes is computed using m o (t)=((E+) 
+ (E-))/2. This mean is subtracted from the original signal 
to get P 1 (t)=P o (t)-m o (t). The procedure is repeated until the 
mean of the envelopes is close enough to zero. If the proce¬ 
dure was repeated n times until to reach the zero-criterion, 
then P m (t) would be first intrinsic mode function (IMF1). 
After the initial intrinsic mode function is found, it is sub¬ 
tracted from the original time series, P 0 (t), and the procedure 
is repeated to locate the second intrinsic mode function. The 
above process is repeated until satisfies above two condi¬ 
tions. The first intrinsic mode function corresponds to the 
highest frequency component of the signal and lower fre¬ 
quency components are extracted in the subsequent intrinsic 
mode functions. The last intrinsic mode function always rep¬ 
resents the climate average, which remains almost constant. 

y(t) = 2" =1 Cj(t) + r n (t) (1) 

here, n is the number of empirical modes, r n is the final resi¬ 
due and c j (t) is the decomposer component. 

Ensemble Empirical Mode Decomposition 
(EEMD) technique 

A noise-assisted data analysis method is Ensemble Empirical 
mode Decomposition (EEMD), represents a major improve¬ 
ment of the Empirical mode Decomposition method, elimi¬ 
nating largely the mode mixing problem and preserving the 
physical uniqueness of decomposition (Wu et al., 2009). The 
principle of Ensemble Empirical mode Decomposition is to 


add white noise, which populates the whole time-frequency 
space uniformly with the constituent components of different 
scales separated by filter bank (Flandrin et al., 2004; Wu and 
Huang, 2004). The Ensemble Empirical mode Decomposi¬ 
tion process is explained as follows: 

1. Add a white noise series to the targeted data set. 

2. Decompose the data with added white noise into in¬ 
trinsic mode functions using Empirical mode Decom¬ 
position. 

3. Repeat step 1 and step 2 again and again, but with dif¬ 
ferent white noise series at each time. 

4. Obtain the (ensemble) means of corresponding intrin¬ 
sic mode functions of the decompositions as the final 
result. 

The number of trials in the ensemble N, has to be large. In 
this study, alpha was set to 0.21 and N was set to 210. 

c,(t) = lim N ^ ^2k=i{cj,k(t) + ar k (t)} 

CO 

In which 

Cj^Ct) + ar k (t) (3) 

Is the k th realization of the j th IMF in the noise added signal, 
a is the standard deviation of the added noise, r k (t[) is the 
residual after extracting first k intrinsic mode function com¬ 
ponents. Here, we set a value is 0.19 and ensemble (N) value 
is 200. 


RESULTS 

Figure. 1 shows that the mean precipitation over India is 
completely decomposed into thirteen intrinsic mode func¬ 
tions can be extracted using empirical mode decomposition 
(EMD) technique, but only eleven intrinsic mode functions 
are shown to characterize the most important components. It 
can be easily observed that all intrinsic mode functions ex¬ 
hibit slow varying amplitudes and frequencies. Each intrinsic 
mode function component denotes the variation of different 
timescales. In order to investigate the gross characteristics of 
oscillations with dominant periods, we applied Lomb Scargle 
(L-S) periodogram analysis to the precipitation data; each in¬ 
trinsic mode function and resultant amplitude spectral plots 
are plotted at the right side of the Figure 1. Compared to 
the other methods, the Lomb Scargle (L-S) method weights 
the data on a per point basis instead of a per time interval 
basis (Press et al., 1992). This technique is equivalent to a 
pure harmonic least-square analysis. The advantage of this 
method is that the input data do not have to be evenly spaced 
in (Scargle et al., 1982; Press et al., 1992). Semi-annual and 
annual oscillations are observed in the second mode of in¬ 
trinsic mode function (IMF2), and third mode of intrinsic 
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mode function (IMF3), respectively. These two oscillations 
with maximum amplitudes of periodicities about 6 and 12 
months are observed. The fifth mode intrinsic mode function 
(IMF5) corresponds to quasi-biennial oscillation (QBO), its 
periodicity between 20-34 months and the maximum ampli¬ 
tude at around 24 months. The seventh intrinsic mode func¬ 
tion (IMF7) has broad periods (3.8 to 6 years) and it cor¬ 
responds to the El-Nino southern oscillation (ENSO) cycle. 
The IMF7 dominates with maximum amplitudes at ~4years 
and ~6 year periods. The tenth mode of intrinsic mode func¬ 
tion (IMF 10), a clear peak between 9 and 11 years, and the 
maximum peak at around 10-year exists, and it corresponds 
to the solar cycle. The eleventh mode intrinsic mode function 
(IMF11) corresponds to the inter-decadal oscillations (IDO) 
and the period oscillates between 18-22 years periods. All 
these oscillations are observed with 90% confidence level. 
Agnihotri et al. (2011) also have reported this inter-decadal 
(16-30 years) variability Total Solar Irradiance (TSI) and 
Indian rainfall datasets. Except for 6-month and 12-month 
oscillations, the remaining oscillations quasi-biennial oscil¬ 
lation (QBO), El-Nino southern oscillation (ENSO), solar 
cycle, and Inter-decadal oscillation are not clearly seen in 
the original time series data set. 

From Figure 1, frequent occurrence of mode mixing, which 
is defined as a single intrinsic mode function either consist 
of widely desperate scales, or a signal of a similar scale re¬ 
siding in different intrinsic mode function components. To 
overcome the scale separation issue without introducing a 
subjective intermittence test, a new noise-assisted data anal¬ 
ysis (NADA) method, known as ensemble empirical mode 
decomposition (EEMD), which defines the true IMF com¬ 
ponents as the mean value of an ensemble number of trials 
(Wu and Huang et al., 2009). The ensemble empirical mode 
decomposition (EEMD) method detailed procedure given in 
Section 3. 

This method extracts the input precipitation data to 13 in¬ 
trinsic mode functions, but we show only the first 10 intrin¬ 
sic mode functions here in the figure 2 for the clarity of the 
most important components. The semi-annual and annual 
oscillations are found in the modes of IMF2 and IMF4. The 
amplitude of the semi-annual oscillations is smaller than an¬ 
nual oscillations. IMF5 is a mode with a dominant period 
of 24-28 months, but the peak value is at about 26 months, 
and this mode indicates quasi-biennial oscillation (QBO). 
IMF6 has associated with the El-Nino southern oscillation 
(ENSO) cycle about a 4-6 year period. The average period 
of IMF7 and IMF8 corresponds to 132 and 216 months, and 
these represent to solar and inter-decadal oscillations. These 
modes all fall above the confidence interval and therefore are 
significant. It is worth mentioning here that the decomposed 
oscillations are fixed intervals than the empirical mode de¬ 
composition (EMD). This is most likely due to the adding 
the white noise to the Indian meteorological department 


(IMD) of precipitation data. 

Figure 3 and 4 provide a more detailed look at the iterations 
in empirical mode decomposition (EMD) and ensemble em¬ 
pirical mode decomposition (EEMD) methods using precipi¬ 
tation over India. Box plots describe the statistical distribu¬ 
tion of iterations for each intrinsic mode function. The box 
plots identify the five important statistics on two methods. 
The iterations are presented in the form of a box plot or 
box-and-whisker plot showing the minimum, median (q 0 5 ) 
quartile range (box), lower quartile (q 025 ), the upper quartile 
(q 0 75 ), and 5 th and 95 th percentiles (whiskers) across the grid 
cells over the two methods of precipitation dataset. Moreo¬ 
ver, while in the empirical mode decomposition (EMD) case 
the total number of iterations is 740, in the case of ensemble 
empirical mode decomposition (EEMD) is 515 iterations. 
It is clear that the ensemble empirical mode decomposition 
(EEMD) method provides less shifting iterations, less com¬ 
putational time for the given time series dataset. 

We further applied Morlet wavelet analysis for six dominant 
periods of intrinsic mode functions in figure 2 are shown in 
figure 5. In the second intrinsic mode function shows the 
presence of six-month wave, and the maximum amplitude 
at about ~5 mm. In the fourth intrinsic mode function period 
covers annual periods ~12 months. The annual amplitude 
(IMF4) is greater than semi annual oscillations (IMF2), fifth 
intrinsic mode function period covers 23-33 months, and the 
maximum peak at around 28 months is quasi-biennial os¬ 
cillation (QBO), the maximum amplitude of 28 months is 
observed during the periods at 1915-1920, 1960-1972, and 
2003-2009.The sixth intrinsic mode function period covers 

4- 6 year periods is El-Nino Southern Oscillation (ENSO) 
and extends nearly six months period during the observa¬ 
tion period, where the maximum amplitudes are observed 
around 1935-1960, 1995-2005. The dominant oscillation of 
9-14 years is clearly observed in IMF7. The maximum am¬ 
plitude is at about 10.8 years during the period of 1940-1960. 
This intrinsic mode function represents a solar cycle. The 
long period of oscillations 16-23 years is observed in IMF8 
is Inter-decadal oscillation. The amplitudes are observed in 
quasi-biennial oscillation (QBO), El-Nino southern oscilla¬ 
tion (ENSO), and Annual oscillations in the wavelet analysis. 

DISCUSSION 

The empirical mode decomposition and Ensemble empirical 
mode decomposition using to investigate the natural periodic 
oscillations into two sets over India, first one is semi annual, 
annual oscillations is small amplitudes and other oscillations 
are Quasi bi-ennial oscillation and El-Nino southern oscil¬ 
lation is coincided with the IMF1 with an average period of 
2.7 years and second mode IMF2 is a dominant period of 

5- 6 years and the long periods are solar cycle, Inter decadal 

■ 
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oscillation are also observed in IMF3 is associated with sun¬ 
spot cycle of about 11 years and fourth IMF is about 20-24 
years in all India rainfall using EMD method by Iyenger et 
ah, (2005). The importance of these oscillations are changed 
the variability of rainfall. 

CONCLUSIONS 

In the present Study, we have investigated the natural pe¬ 
riodic oscillations using India Meteorological Department 
(IMD) precipitation data during the period from 1901-2010. 
We examine the semi annual oscillation (SAO), annual os¬ 
cillation (AO), quasi bi-ennial oscillation (QBO), El~Nino 
southern oscillation (ENSO), inter-decadal oscillation (IDO) 
periods using Empirical mode decomposition (EMD) and 
Ensemble Empirical mode decomposition (EEMD) tech¬ 
nique. We clearly found mixed mode oscillations in Em¬ 
pirical mode decomposition. The Empirical mode decom¬ 
position method shows a strong mixture of modes, where 
Intrinsic mode function dramatically changes with no valid 
reason. This instability can have a dramatic effect on the fur¬ 
ther study of any signal. According to the obtained results, 
the Ensemble Empirical mode decomposition improves the 
precipitation intrinsic mode functions, and offers a simple 
approach for the stable prediction of non-stationary data. As 
a future work, it would explore the possibility of employing 
different aggregation methods as well as performing an ad¬ 
dition and more significant test that exposes more reliable 
results, may be considering for other datasets. 
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Figure 1: Time series of monthly mean India Meteorological 
Department (IMD) precipitation during 1901 to 2010. Intrinsic 
mode function (IMF) components extracted from Empirical 
mode decomposition (EMD) method, the first to eleven in¬ 
trinsic mode functions are shown in left column. Correspond¬ 
ing Lomb-Scargle periodograms are shown in right column. 
Dashed horizontal line indicates 95% confidence level. 
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Figure 2: Time series of monthly mean India Meteorological 
Department (IMD) precipitation during 1901 to 2010. Intrinsic 
mode function (IMF) components extracted from Ensemble 
Empirical mode decomposition (EEMD) method, the first to ten 
intrinsic mode functions are shown in left column. Correspond¬ 
ing Lomb-Scargle periodograms are shown in right column. 
Dashed horizontal line indicates 95% confidence level. 
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Jyothi et.al.: Natural periodic oscillations extracted in the precipitation using empirical mode decomposition and ensemble empirical ... 



Figure 3: Box plot of Empirical mode decomposition (EMD) 
performance measures for the monthly precipitation dataset. 
The vertical axis is number of iterations. 




Figure 5: Contours of wavelet intensities of frequency with 
time for the precipitation datasets and its six intrinsic mode 
functions (IMFs) are (2, 4, 5, 6, 7, and 8) as shown in Figure.2. 


Figure 4: Box plot of Ensemble Empirical mode decomposi¬ 
tion (EEMD) performance measures for the monthly precipita¬ 
tion dataset. The vertical axis is number of iterations. 
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